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We report the results of a search for a charmonium-like state produced in the process 77 — > ojj/ifj 
in the 3.9-4.2 GeV/c 2 mass region. We observe a significant enhancement, which is well-described 
by a resonant shape with mass M = (3915±3±2) MeV/c 2 and total width Y = (17 ± 10 ±3) MeV. 
This enhancement may be related to one or more of the three charmonium-like states so far reported 
in the 3.90-3.95 GeV/c 2 mass region. 

PACS numbers: 13.25.Gv, 13.66. Be, 14.40.Pq, 14.40.Rt 



Many new charmonium-like states have been discov- 
ered by the B-factory experiments, typically via a promi- 
nent hadronic decay to a known charmonium state, such 
as J/ip, ip(2S) or Xd- Some have attracted particular 
interest because of their net electric charge Three 
of the new (neutral) states were discovered by Belle in 
the 3.90-3.95 GeV/c 2 mass region. The X(3940) was 
found in the e + e~ — > J jvbX double charmonium produc- 
tion process, with a prominent decay to the DD* final 
state 0. The F(3940) was observed in the B decay pro- 
cess B~ -> r(3940)if- with F(3940) -> uJ/ij} UH, 
and is a candidate for an exotic state, such as a hybrid 
meson (ccg), or a D*D* bound state Jl. The Z(3930) 
was found in the 77 — > DD process [6|, and is usually 
identified with the x C 2(2P). These three states appear 
in different production and decay processes, and are usu- 
ally considered to be distinct particles, however there is 
no decisive evidence for this. The interpretation of these 
states has been discussed by many authors: see, e.g., 
Ref. 0. 

It is important to search for a signature of the F(3940) 
or any other resonant state contributing to two-photon 
production of uiJ/ip. This final state is the lightest com- 
bination of two vector mesons with definite C-even and 
1 = quantum numbers that can be produced in two- 
photon processes via a hidden-charm state. In this paper 
we present measurements of the 77 — > ujJ/ip process in 
the 3.9-4.2 GeV/c 2 mass region, in which we observe 
a resonant enhancement. The signal is from the two- 
photon process e + e~ — > e + e~ujj/ij) in the "zero-tag" 
mode, where neither the final-state electron nor positron 
recoiling from photon emission arc detected. 

We use experimental data recorded with the Belle de- 
tector Q at the KEKB e + e~ asymmetric-energy (3.5 on 
8 GeV) collider @, corresponding to an integrated lu- 
minosity of 694 fb _1 . The data are accumulated mainly 



on the T(45) resonance (y/s = 10.58 GeV) and 60 MeV 
below it. A small fraction of data from different beam en- 
ergies near 10.36 GeV (the T(35) mass) and 10.87 GeV 
(the T(55 l ) mass) is also included in the sample. 

A comprehensive description of the Belle detector is 
given elsewhere §|. Charged tracks are reconstructed 
in a central drift chamber (CDC) located in a uniform 
1.5 T solenoidal magnetic field. The z axis of the de- 
tector and the solenoid is along the positron beam, with 
the positrons moving in the —z direction. Track trajec- 
tory coordinates near the collision point are measured 
by a silicon vertex detector (SVD). Photon detection 
and energy measurements are provided by a CsI(Tl) 
electromagnetic calorimeter (ECL). A combination of 
silica-aerogel Chercnkov counters (ACC), a time-of- flight 
counter (TOF) system consisting of a barrel of 128 plas- 
tic scintillation counters, and specific ionization (dE/dx) 
measurements in the CDC provides K/tt separation for 
charged tracks over a wide momentum range. The mag- 
net return iron is instrumented to form a Kr, detec- 
tion and muon identification (KLM) system that detects 
muon tracks. 

Signal candidates arc triggered by a variety of track 
triggers that require two or more CDC tracks associated 
with TOF hits, ECL clusters, a total energy deposit in 
the ECL above a threshold (0.5 GeV), or a muon track 
in the KLM detector. In addition, events with a total 
ECL energy above 1.1 GeV are triggered by a separate 
logic. Because of the presence of a lepton pair in the 
final state of the signal processes, a combination of the 
above triggers provides a high overall trigger efficiency, 
(98 ± 2)%. 

We select signal event candidates by reconstructing 
all the final state particles from u — > ir + 7r~ir° and 
J I if) — > l + l~ (I = e or fi). Twelve selection criteria are 
imposed: (1) there are just 4 charged tracks with trans- 
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verse momentum p t > 0.1 GeV/c originating in the beam 
collision region; (2) the net charge of the tracks is zero; 
(3) none of the tracks is identified as a kaon (we require 
a likelihood ratio £{K)/(£(K) + £(tt)) < 0.8, which is 
satisfied by 99.6% of pions but only 5% of kaons, for mo- 
menta below 0.8 GeV/c); (4) there is a net-charge-zero 
combination of two tracks whose invariant mass is in the 
J/tp mass region, |M(2 tracks) - Mj/^\ < 0.2 GeV/c 2 , 
where Mj/^, = 3.0969 GeV/c 2 and we assume the pion 
mass for each of the two tracks; (5) there is one or more 
neutral pion candidate formed by a mass-constrained fit 
of two photons, with \ 2 < 4; (6) the number of 7r°'s 
with p t > 0.1 GeV/c must not exceed one in an event. 
If there is a ir° satisfying the pt condition, it is accepted 
as the only ir° candidate in that event. If no 7r° satisfies 
the pt condition, we retain all the ir° candidates (below 
Pt < 0.1 GeV/c) at this stage; (7) events in the kinematic 
region of initial-state-radiation processes, e + e~ — > ■yX, 
where the photon is emitted very close to the direction 
of the incident e~ beam, are eliminated. We reject events 
satisfying the following condition on the z-component of 
the laboratory momentum for the final-state system X 7 
P z < (Ml - 49.0 GeV 2 /c 4 )/14.0 GcV/c 3 + 0.6 GeV/c, 
where P z and M5 are the z-component of the momen- 
tum and the invariant mass of the system constructed 
from the four tracks and a neutral pion candidate, re- 
spectively. 

For further analysis we examine only events with 
W < 4.3 GeV, where W is defined by W = M 5 - 
M(l + l~) + Mj/^, using a refined two-lepton invariant 
mass (M based on the lepton flavor identified by 
the following criteria: (8) if either of the tracks is identi- 
fied as an electron, based on the ECL energy deposit, the 
tracks arc identified as e + e~. Otherwise, if cither track 
is identified as a muon based on KLM information, the 
tracks arc identified as p + fi~ . An event that fails both 
tests is rejected. If one or more photons with energy be- 
tween 20 and 200 MeV are found within 3° of either the 
e + or e~ track, the energy of the most energetic photon 
near the track is added to the track momentum. 

Following this correction, (9) we refine the J/tp 
selection with a more stringent requirement for the 
lepton-pair invariant mass, 3.07 GeV/c 2 < M(l + l~) < 
3.12 GcV/c 2 ; (10) we suppress ifj{2S)n° events with the 
mass difference requirement, \M(l + l~ir + ir~) — M — 
0.589 GeV/c 2 | > 0.01 GcV/c 2 ; (11) to select uj can- 
didate, a condition on the ir + n~ir° invariant mass, 
0.753 GeV/c 2 < M(3vr) < 0.813 GcV/c 2 , is imposed. 
If there are multiple uj candidates due to multiple 7r°'s 
in an event, we choose the one with the smallest \ 2 m 
the 7T° mass constrained fit, in order to avoid multiple 
entries in the final ujj/ip mass spectrum. Finally, (12) 
we require transverse momentum balance for the 5-body 
system, < 0-1 GeV/c, where p\ is the momen- 

tum of a particle in the e + e~ cm. frame, in the plane 
perpendicular to the beam direction (loj . 
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FIG. 1: (a) The M(l + l~) distribution just after the dilepton 
selection, (b) The M (l + 1~ tt + tt~ ) - M(l + l~) distribution af- 
ter the tight J/tp selection. Events between the arrows are 
rejected as consistent with ip{2S) production. 

Figures 1(a) and 1(b) show the distributions of 
M(l + l~) just after requirement (8) and the mass dif- 
ference M(l + l~ir + ir~) — AI(l + l~) just after requirement 
(9), respectively. The ip{2S) contribution is effectively 
removed by criterion (10) above. 

The main background process is multi-pion produc- 
tion from two-photon processes. However, after all se- 
lection requirements are applied, non-LoJ/ip backgrounds 
are rather small, as shown in the scatter plot in Fig. 2(a) 
for the samples where all the selections except those for 
M(l + l~) and M(3tt) arc applied. Figures 2(b) and 2(c) 
show the distributions of M(l + l~) and M(37r), respec- 
tively, in the selection bands for the opposite-side par- 
ticle; for clarity, we exclude events somewhat below the 
uJ/ip threshold, W < 3.85 GeV. 

In Figs. 2(b) and 2(c), the experimental M(l + l~) and 
M(3ir) distributions are compared with those from the 
signal Monte Carlo (MC) events, which are generated 
assuming spin-parity (J p ) and mass (W) of the ujJ/tp 
system to be + and 3.93 GeV/c 2 , respectively. Details 
of the signal MC generation are given below. We confirm 
that the experimental mass distributions are consistent 
with those of signal MC events. 

We find that there are two events in the signal region 
with multiple to candidates, out of 73 events in total; 
we choose only one combination in each event, according 
to criterion (11). The fraction is consistent with the 1- 
2% multiple candidate rate expected from the signal MC 
sample. 

We show the W distribution for the final 77 — ^ 
ojJ/'ijj candidate events in Fig. 3. There is a prominent 
resonance-like peak around 3.92 GeV. It is far above the 
non-uiJ/t/j background contribution, which is estimated 
from the events in the oj and J/tp mass sidebands (shown 
as shaded histograms for comparison); we define eight 
sideband regions in the plane of Fig. 2(a) with the same 
dimensions as the signal region, i.e., each region centered 
at 3.035 GeV, 3.095 GeV or 3.155 GeV with a width 
of 0.05 GeV in the M(l + l~) direction and centered at 
0.693 GeV, 0.783 GeV and 0.873 GeV with a width of 
0.06 GeV in the M(37r) direction, and average the dis- 
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FIG. 2: (a) Scatter plot of M(3tt) vs. M(l + l~) for the exper- 
imental data, with all the other cuts applied. The rectangle 
shows the signal region, (b) The M(l + l~) and (c) M(3tt) 
distributions with all the other cuts applied, and requiring 
W > 3.85 GeV. Entries due to multiple cj candidates in a 
given event are included. Points with error bars show the 
data; the dashed histograms show signal MC events, normal- 
ized to the data yield in the selection area, with the scaled 
sideband yield subtracted. The arrows show the selection re- 
gions. 

tribution over the eight regions. We modify the W value 
of each sideband event plotted in Fig. 3, shifting it by 
the difference between the sum of mass coordinates of 
the central point of the signal region (3.878 GeV) from 
that of the sideband region where the event is found, for 
comparison to the signal-event distribution. 

Figure 4(a) shows a scatter plot of the transverse mo- 
mentum balance vs. W after requirement (11). A promi- 
nent concentration of events near W = 3.89 — 3.95 GeV 
and | YlPt I < 0.05 GeV/c is visible; a comparison of the 
| J2p*\ projection with signal MC is shown in Fig. 4(b). 
Based on these results, and the shape in W (Fig. 3), we 
conclude that the concentration of events is due to a res- 
onance formed in two-photon collisions. 

The W distribution for the final candidate events 
is fitted by an incoherent sum of resonant and back- 
ground components. We adopt an S-wave Breit-Wigner 
function with a variable width for the resonant com- 
ponent, (2N r /tt)M 2 T'/{(W 2 - M 2 ) 2 + M 2 T' 2 )} and 
r' = T(p*/pq), where p* is the momentum of the two- 
body decay to ujj/ip, in the rest frame of a parent par- 
ticle of mass W] p* is the value for W = M The 
nominal mass (M), width (r) and yield parameter {Nr) 
are treated as fit parameters. 

We represent the background component by a 
quadratic function of p* that vanishes at the nominal 
uJ/ip threshold, M th = 3.8796 GeV/c 2 . We also add a 
constant term, to represent the high W tail, which, as 
the sideband study suggests, is dominated by xion-uJ/tp 
events. The sum of the two components has a functional 
form, {{ap* + bp* 2 ) + c}9(W - Mth), where 0(x) is a unit 
step function that is non-zero only for x > 0. The param- 



eters a, b and c are floated within the constraint that each 
of the two background components must be non- negative 
throughout the fitting region. 

The fit takes into account the W resolution in the mea- 
surement, which is approximated by a double-Gaussian 
function from the signal MC events (59% of the sig- 
nal has a resolution a of 4.5 McV, while the remainder 
has a = 16 McV with the peak position displaced by 
—4 MeV). We perform an unbinned maximum likelihood 
fit in the region 3.875 GeV< W < 4.2 GeV. The signal 
candidates with the smallest W are the two events with 
W between 3.879 and 3.880 GeV. 

The W dependences of the efficiency and luminosity 
function are taken into account in the fitting function. 
The efficiency is determined using signal MC events as 
described in detail later. We use the W dependence of 
the efficiency for J p = + for the nominal fit. Between 
the threshold and 3.96 GeV, the VT-dependence is weak: 
the efficiency varies by 10% only, and has a minimum 
near W = 3.92 GeV. 

The obtained resonance parameters for the mass and 
the width are as follows: 

M = (3915±3±2) MeV/c 2 , 
T = (17 ±10 ±3) MeV, 

where the first and second errors are statistical and sys- 
tematic, respectively. The estimated yield from the res- 
onant component in the fit is 49 ± 14 ± 4 events in the 
region below 4.2 GeV. The statistical significance of the 
resonant peak is 7.7cr, which is determined from the dif- 
ference of the logarithmic likelihoods, — 21n(L/Lo), tak- 
ing the difference of the number of degrees of freedom in 
the fits into account, where L$ and L are the likelihoods 
of the fits with and without a resonant component, re- 
spectively. The relevant fit curves are shown in Fig. 3. 
The x 2 °f the nominal fit, determined using 10-MeV- 
width binning in the range 3.85-4.2 GeV, is 27.3, for 29 
degrees of freedom. 

The systematic errors quoted above are determined 
from a study of alternate fits: we use a Breit-Wigner 
function with a constant width; we enlarge the invariant- 
mass resolution by 20% (an over-estimate of the data- 
Monte Carlo difference allowed by the fit); we change 
the upper limit of the fit region in W to 4.1 GeV and 
4.3 GeV, respectively. The changes in the central values 
of the corresponding resonance parameter are combined 
in quadrature. We also take into account the uncertainty 
of the mass scale, estimated to be 1 MeV/c 2 , in the mea- 
surement of M. There is no significant change in the 
parameters if J p = 2 + is assumed; the changes of mass 
and width are less than 0.1 MeV/c 2 and 0.3 MeV, re- 
spectively. The resonant contribution for the J p = 2 + 
assumption is 1.0 event smaller than that for J p = + . 

The efficiency for selecting 77 — ¥ loJ/i/j events is de- 
termined using signal MC events generated by TREPS 
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FIG. 3: The W distribution of the final candidate events (dots 
with error bars). The shaded histogram is the distribution of 
non-tj J / 'tp backgrounds estimated from the sideband distribu- 
tions. The bold solid, thinner solid and dashed curves are the 
total, resonance and background contributions, respectively, 
from the standard fit (see the text). The dot-dashed curve is 
the fit without a resonance. 
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FIG. 4: (a) Scatter plot of pt balance vs. W for the final 
candidate events in which only requirement (12) is omitted, 
(b) The projection onto the pt balance axis for events with 
W < 3.95 GeV. The dashed histogram is the expectation 
from signal MC events, normalized to the number of signal 
candidates in the selected region. The pt balance requirement 
is indicated by the arrow. 



code [ll|. We generate 10 5 MC events for both e + e~ and 
decays of J/ip, at nine different W points between 
3.89 and 4.15 GeV. The efficiency for the signal process 
at W = 3.92 GeV is determined to be (1.85 ± 0.20)% 
((1.26 ± 0.14)%) for the J p = 0+ (2+) assumption; the 
efficiency is defined for the range Q 2 < 1.0 GeV 2 , for each 
incident photon. We assume production in helicity-2 for 
J p = 2 + [6] and decay to ujj/ijj in an S-wave for both 
+ and 2 + . The other two possible J p assumptions, 0~ 
and 2 _ , are similar and give an efficiency close to that 
for J p = 0+. 

Based on the efficiencies calculated for the two J/ip 
decay modes, the fraction of signal in the e + e~ mode is 
expected to be 36%. This is consistent with the fraction 
in the data: 27 J/ip — > e + e~ events among the 73 signal 
candidates. 

Sources of systematic errors in the efficiency determi- 
nation and their contributions are listed in Table I. We 
confirm that the inefficiency due to each of the particle 
identification cuts, (3) and (8), is very small, less than 



TABLE I: Sources and sizes of systematic error in the effi- 
ciency determination 



Source 


Syst. error (%) 


Trigger efficiency 


2 


Track reconstruction 


4 


tt® reconstruction 


3 


Particle identification cuts 


2 


Effect of background hits 


3 


J/ip selection 


3 


uu selection 


6 


W-dependence, effect of background, etc. 


3 


Luminosity function, integrated luminosity 


5 


Total 


11% 



1%, for signal events. The uncertainties in the efficiencies 
of the invariant mass cuts are estimated by varying the se- 
lection regions near M(l + l~) = Mj/^ and M(37r) = M u 
by ±20% in the MC. We sum the uncertainties in quadra- 
ture, and find 11% in total. 

Treating the observed structure as a resonance denoted 
by A"(3915), we derive the product of the two-photon 
decay width and the branching fraction to wJ/ip, using 
the yield parameter Nr from the fit and the selection 
efficiency. We obtain 

r 77 (V(3915))i3(V(3915) uJ/ip) 

f (61 ± 17 ± 8) cV for J p = 0+ 
~ \ (18 ± 5 ± 2) eV for J p = 2+ helicity-2. 

Based on this result, and the measured width T, the 
product of the two partial widths of the A"(3915), 
T 11 (X)T UJ j^(X) is of order 10 3 keV 2 . If we assume 
j 0(1 keV), typical for an excited charmonium 
this implies T^j/^ ~ 0(1 MeV): a rather large 
value for a charmonium-transition partial width of such 
a state. This value of the product of the partial decay 
widths is roughly compatible with the prediction assum- 
ing the D*D* bound-state model 0]. 

To conclude, we have observed a resonance-like en- 
hancement in the 77 — > uiJ/tp process with a statis- 
tical significance of 7.7er, which contains 49 ± 14 ± 4 
events in the peak component. The mass and width have 
been measured to be M = (3915 ± 3 ± 2) McV/c 2 and 
T = (17 ± 10 ± 3) MeV, respectively These values are 
consistent with those of the V(3940), which is seen in the 
ujj/ip final state 0,0]; an d close to those of the Z(3930), 
which is seen in 77 — > DD Q . 
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